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ABSTRACT 



A forecast model proposed by 'Arnason as being capable of producing 
realistic forecasts of displacement of systems at and below 500 mb is 
investigated. The prediction equation. 



is applied to test cases at 1000, 850, 700 and 500 mb with various values 
of the parameter ki* The effect of the value of the parameter on the 
amount of divergence present in the model is discussed. Root -mean-square 
errors for 2U-hour and US-hour forecasts with various values of at 
each level are presented. Optimum values of k^^ for the cases tested are 
selected as 3.5> 2.75 and 0.9 for the levels 850, 700 and 500 mb respec- 
tively. It is found that with proper values for k^ the model will provide 
good forecasts of system movement at each level tested, but as expected, 
changes in the intensity of systems are not correctly forecast. In all 
cases values of root-mean- square error of forecasts are less than those 
obtained by applying the "Helmholtz" barotropic model. A geostrophic 
version of the model is also tested and is found to produce results simi- 
lar to those obtained for the stream-function version except that the 
geostrophic "blow-up" of low-latitude high pressure cells is observed. 

The writer wishes to express his appreciation for the guidance and 
encouragement given him by Professor George J. Haltiner of the U. S. 

Naval Postgraduate School in the course of this investigation. In addi- 
tion, appreciation is expressed to Mr. Geirmundur 'Arnason, on whose work 
this paper is based, and to personnel of the Fleet Numerical Weather 
Facility for many instances of assistance. 
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Introduction 



Since the initation of operational numerical weather forecasts by 
the Joint Numerical Weather Prediction Unit in 1955^ the only forecast 
model which has stood the test of daily operational use has been the 
barotropic model first proposed by Rossby in 1939. The barotropic model 
remained in its original form until 1958 when Wolff [1] discovered large 
systematic errors due to improper forecast of very long waves. After 
these errors had been empirically corrected by Wolff, it was discovered 
by Cressman [.2] that these long-wave errors could be reduced by the in- 
troduction of divergence into the barotropic model. This resulted in 
the divergent or "Helmholtz” version of the barotropic model which is in 
general use at this time. 

With the improved barotropic model, high quality forecasts have been 
available for the 500-mb level. This barotropic model, however, is by 
nature restricted in its use to the 500-mb level and, in addition, is ca- 
pable mainly of forecasting the movement of systems, and not in general 
changes in intensity. Occasionally attempts have been made to replace 
the barotropic model in operational use with some type of baroclinic 
model which would obviate these shortcomings. These attempts have so far 
met with failure largely due to the inability to control the forecast of 
cyclogenesis. 

’Arnason has suggested that, in view of the great utility of the 
500-mb forecasts produced by the barotropic model, it might be worth- 
while temporarily to set aside the problem of predicting cyclogenesis 
and to try to develop a model which would possess the same skill of fore- 
casting movement at other levels which the barotropic has at 500 mb. To 
this end, *Arnason has proposed a model [ 3 ] which maintains many of the 
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characteristics of the barotropic, but which is more adaptable to other 
levels. 

This model is based on a single-layer compressible fluid with a 
stratified density distribution, bounded by two constant-pressure sur- 
faces. The theoretical development results in an expression for the 
horizontal wind divergence 

J.vW . tV- V^'l (1) 

J . 

Introducing this expression into the vorticity equation, 

IS + V- V7-7 c)>v V = O , (2) 

gives the prediction equation for the "stratified” model 

:>!. o. ( 3 ) 

,u ' / 

In order to apply this equation to forecasting, 'Arnason has suggested, 
as one alternative, the use of a stream function for computation of wind 
and vorticity and the elimination of ^ through a geostrophic assump- 
tion. This gives 

^\V-V7 - £7. V ' V ? - 0, (h) 

which is of the same form as the prediction equation for the divergent 
or "Helmholtz" barotropic model, but with an additional divergence term 
involving V and with the substitution of the stability-dependent term 
g* for g. 

Although the forms of the prediction equations are similar, there 
are several basic differences between the "stratified" and other barotrop- 
ic models. By the method of perturbation analysis 'Arnason has obtained 
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the meteorologically significant root of the frequency equation for the 
stratified model 



C = 




( 5 ) 



Although this equation is again similar to the corresponding equation 
for the Helmholtz barotropic model, an important difference arises in 
the fact that the basic current V does not have a multiplying factor. 
By empirically varying the value of g'H, the magnitude of the ^ term 
may be changed without reducing the effect of the basic current, thus 
allowing control over the speed of movement of systems. As can be seen 
from equation (l) this is equivalent to a variation in the amount of 
divergence present in the model. This provides a method of ad;justing 
the model to yield proper system movement at a number of levels. 

The stratified model as proposed by 'Arnason seems to provide a 
promising approach to numerical forecasting at and below the 500-mb 
level. It is the purpose of this paper to apply the model to a number 
of test cases and determine its effectiveness. 
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. Application of the Model to Forecasting 

In the application of the prediction equation ( 3 ) to actual fore- 
casts there are several choices to be made. The parameter V has not 
been explicitly defined, but may be taken to be any kind of space- 
averaged or zonal wind. Also the value of g*H is left to be determined 
empirically. 

By examination of the prediction equation ( 3 )j it can be seen 
that the additional "divergence" term of the stratified model which in- 
volves V essentially represents the horizontal advection of height 
by the mean wind. This immediately suggests the notion of a steering 
current, which may be thought of as a wind at some level or as a mean 
wind in some layer of the atmosphere. Also, the magnitude of the di- 
vergence term is dependent on the difference between the mean wind and 
the wind at the forecast level. Thus it appears reasonable to use as 
a mean wind something close to the 500-mb wind, since this would mini- 
mize the divergence term at SOO mb and allow it to increase with dis- 
tance from 500 mb. This is generally consistent both with observed 
patterns of divergence in the atmosphere and with the need for increas- 
ing the value of the divergence term at lower levels to supplement the 
movement due to vorticity advection alone. 

A second parameter which must be empirically interpreted is the 
factor g'H. The value of g' is a measure of the static stability of 
the basic flow and H is the thickness of the layer considered in the 
model. Since neither is uniquely defined, the combination may be taken 
as any reasonable value which produces proper movement of the systems 
being forecast. 

In order to make a forecast, either or must be eliminated 

h 



from equation (3). 'Arnason has used a geostrophic relation to eliminate 



tion, height fields must first be "balanced" to obtain the stream function 
and then, at the conclusion of the forecast, inverted to obtain the fore- 
cast height field. Another approach would be to compute the vorticity 
and winds directly from the height values using the geostrophic rela- 
tionship, obtaining the forecast equations 



^ , obtaining equation (I 4 ) which involves f alone. To apply this equa- 






n t 




l_2L V ■ vz = o 



(6) 






TT S' 
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3. Procedure 

As mentioned above, in order to apply the stratified model to actual 
forecasting, values must be chosen for the parameters V and g'H. For 
the purposes of the tests presented here, it was decided to use for \/ 
a mean wind computed from a height field obtained by first heavily 
smoothing the 500-mb and 700-mb heights, then computing an intermediate 
field by the equation 

^ .IS .25 i> . (7) 

500 TOO 

The wind obtained by this method is thus a space-mean wind for some 
intermediate level between 500 mb and 700 mb. For the first six one- 
hour time steps, a V computed from the initial data was used; there- 
after a new V was computed from forecast fields every sixth time step. 
In the choice of the parameter g'H, consideration was given to the 
similarity between the 'forecast equation (Ii) and the forecast equation 
for the "Helmholtz" barotropic 




Since it was known that forecasts made with this model at 500 mb are of 
consistently good quality, it was decided to choose g'H such that the 
Helmlioltz term for the stratified model would be the same as the Helm- 
holtz term of the barotropic. Thus in the Helmholtz term, g'H was 
chosen to be 

'H = cr H.P (9) 

where k is a "tuning constant" with a value of four, and H^ is the stand- 
ard height of the 500-mb surface. 

In order that the model be used at levels other than 500 mb, it was 
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decided to include a second "tuning constant" in the final divergence 
term. By varying this second constant the proper amount of divergence 
may be introduced at each level to produce the best forecast. With 
these values introduced, the forecast equation for the stream function 
case becomes 

{ +\V-V'7 V'Vf^O (10) 

^ IHT 

In the same manner, starting with equation (6) a similar equation was 
obtained for the geostrophic version of the model: 

/ ^ 2: +LV'V7 - V- V2: = o (11) 

^ gHTi JT ^ 

The model, in the form of finite difference equations obtained from 
equations (lO) and (ll), was run to li8-hour forecasts for 850 mb, 700 mb 
and 500 mb for several randomly selected days. Various values of were 
used at each level for both the geostrophic and stream-function cases. 

In addition, for the 6 May 1962 case runs were made on 1000-mb data in 
order to determine the feasibility of using the model at that level. 
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U. Computational Methods 

The test computations of the stratified model described in this 
paper were made on a Control Data Corporation l60U computer using a 
program constructed of meteorological subroutines prepared by the Fleet 
Numerical Weather Facility, Monterey. The grid used in the numerical 
solution of the prediction equation was the standard octagonal grid used 
by both the Joint Numerical Weather Prediction Unit at Suitland, Md. and 
the Fleet Numerical Weather Facility. 

Time differencing was accomplished in one-hour time steps using the 
central difference formula 

. ( 12 ) 

where is the computed one-hour change in geopotential. At the 

initial time, T= O , and all integral multiples of twelve, the forward 
difference formula 

<i>r„ = K+ 

was used. 

In cases where a stream function wind was used, the stream values 
were obtained from initial height data through the solution of the bal- 
ance equation as described by 'Arnason [1^1. All initial height fields 
and fields required for forecast verification were obtained from files 
of operational analyses at the Fleet Numerical Weather Facility, 
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5. Results 

Verification scores for each of the cases tested were determined 
by a FNWF program using the equations 



Here A represents the forecast value, B, the observed value and N, the 
total number of grid points, 1977. 

As can be seen from the equations, the pillow is a measure of the 
error in the forecast of the mean height while the RMSE (root mean square 
error) is a measure of the deviation of errors around the mean error. 

Figures 1, 2 and 3 show plots of RMSE as a function of for the 
various cases tested. As can be seen, the value of necessary to pro- 
duce the best forecast (minimum RMSE) varies somewhat from case to case 
at each level. However, the curves are relatively flat in the vicinity 
of the minima so that by selecting a reasonable value of k^^, the quali- 
ty of forecasts from one case to another will vary only slightly from 
the optimum. Table 1 shows mean values of RMSE for all cases obtained 
with various values of k^ at each level. A plot of these values is 
shown in figure U. From this plot of mean error versus k^, values of 
k^^ can be selected which will give the minimum total error. Table 2 
gives a tabulation of the error which would occur in each case tested if 
kj at each level were selected as the value which produces the smallest 
mean error. 

Forecasts for 8^0 mb made from initial data for 20 January 1962 are 



Pillow 




ilk) 



N 




( 15 ) 
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Table 1. Mean Verification Scores of Forecasts With Various Values of 



STREAM-FUNCTION FORECASTS 




kl 


RMSE (ft) 


kl 


RMSE (ft) 


1^1 


RMSE (ft) 


kl 


RMSE (ft) 


(Q 


850 


2 


112 


3 


107 


1; 


108 


3 


lilt 


U 




















P 

o 

sz 


700 


1 


133 


2 


12U 


3 


123 


h 


128 


(M 


500 


0 


lii8 


1 

2 


1U7 


1 


Ihl 


2 


l5o 


CO 


850 


2 


17h 


3 


166 


h 


166 


5 


173 
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1 


189 


2 


171; 
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169 
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173 


i CO 
1 ^ 1 
1 
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i 500 


0 


219 


1 

2 


215 


1 


211; 


2 


218 

i 


i 

GEOSTROPHIC FORECASTS 


' 

1 


kl 


RMSE (ft) 


kl 


RMSE (ft) 


kl 


RMSE (ft) 


kl 


RMSE (ft) 


) 

CO 


850 


2 


130 


3 


123 


h 


125 


5 


13U 


u 




















o 

x: 


700 


1 


137 


2 


127 

» 


3 


126 


1; 


131 




500 


0 


163 


1 

2 


161 


1 


162 


2 


169 


CO 


850 


2 


198 


3 


192 


1; 


193 


5 


200 






















o 

XI 


700 


1 


— 


2 


203 


3 


195 


1; 


198 


CO 






















500 


0 


— — — 


1 

2 


271 


1 


269 


2 


277 
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Table 2. Verification Scores For Each Day Using Optimum Value of 





850 mb 
Rl^SS (ft) 


700 mb 
RMSE (ft) 


500 mb 
RMSE (ft) 




ki=3.5 


ki=2.75 


ki=0.9 


U November 1961 
(2li hr) 


111 


116 


13U 


20 January 1962 
( 2 li hr) 


128 


128 


178 


6 May 1962 
(2li hr) 


87 


88 


127 


U November 196 I 
(li8 hr) 


173 


185 


226 


20 January 1962 
(li8 hr) 


198 


185 


250 


6 May 1962 
(kS hr) 


lii3 


lliO 


168 
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shown in figures 8 through 13. Initial charts and charts for the veri- 
fying time of each forecast are shown in figures 5 through 7. Comparing 
forecasts for the various values of at a given level shows that the 
forecast speed of synoptic waves increases as the value of is increased 
(thus decreasing the value of g’H and increasing the magnitude of diver- 
gence). This is in accordance with the results of equation (5). Further, 
it can be seen that there is very little change in intensity of the vari- 
ous systems as k^^ is varied. This demonstrates the point that the addi- 
tional divergence term of the stratified model results in movement rather 
than development . 

As can be seen from equations (lO) and (ll) a choice of k^^^O would 
reduce the forecast to the barotropic. An examination of figure li re- 
veals that the optimum value of k]^ at 500 mb is close to zero, but that 
the stratified model does show a slight improvement over the barotropic 
with values of k^^ between zero and one. This is as expected, since 500 
mb is generally accepted to be at or near the level of minimum divergence, 
and thus the contribution of divergence to the change in the vorticit 3 '’ 
pattern at this level should be small. 

At levels other than 500 mb, where the barotropic model is less ade- 
quate, the values of k^^ for which the minimum RMSE is reached are greater. 
For the 700-mb level optimum values for the cases tested range from k]^=2 
to kx=Uj while at 850 mb they range from 3 to 5. These results can be 
accounted for by the vertical distribution of velocity divergence which 
generally increases with distance from the 500-mb surface. 

By use of equation (1) and the chosen values for g’H given by equa- 
tion (9), an expression for the divergence at a given time may be ob- 
tained: 
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(16) 



- 4-f f A± t k,\v- Vf’ 

% 

In figures lU, 1 $ and l6 divergence charts at 850, 700 and 500 mb as com- 
puted from equation (l6) for initial data of 20 January 1962 are shown. 
These charts show clearly the vertical distribution of divergence, with 
values as large as + 5 x 10"^ sec”^ at 850 mb, + 3 x 10“^ sec ^ at 700 
mb and values very nearly zero at 500 mb. 

An examination of the divergence fields relative to their correspond- 
ing height fields reveals that lines of zero divergence very closely fol- 
low the short-wave trough and ridge lines of the corresponding height 
fields. Between the short-wave trough and ridge lines are located the 
centers of maximum convergence and divergence, with divergence ahead of 
ridges and convergence ahead of troughs, as might be expected. 

For the case of 6 May 1962, 1000-mb forecasts were made. From 
figure 3 it can be seen that the curves of RI'ISE versus k]_ for 1000 mb 
are almost exactly parallel to the curves for 850 mb with an RMSE about 
ten feet greater at 1000 mb than at 850 mb. The fact that the curves 
are nearly parallel suggests that the increased error at 1000 mb is due 
to a greater change in the intensity of systems at this level than at 
850 mb. The movement of established systems which appear on both initial 
charts is nearly the same at both levels j thus which controls the mag- 
nitude of the additional movement term should reasonably be chosen the 
same for the two levels. The verifying charts for the forecasts in 
question, figures 21, 22, 25 and 26, show this to be the case. Forecasts 
at both levels have the major systems placed very close to their actual 
positions at verifying time, but the verifying charts at 1000 mb differ 
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more in pattern due to the variation in intensity change of different 
systems over the forecast period. 

In addition to the forecasts made for the three days using stream 
functions obtained from the initial height fields, for two of the days 
forecasts were made directly from the height fields using equation (11). 
As can be seen from figures 27 through 30, the geostrophic forecasts 
thus obtained showed great similarity to the forecasts made from the 
stream function. However, as expected, the geostrophic version appears 
to suffer from the major error common to geostrophic forecasts - the 
"blow-up" of lovj- latitude high-pressure areas. As can be seen, the i;8- 
hour geostrophic forecast increases the central height of the 500-mb 
high off the east coast of the United States to over 19,500 feet, while “ 
the stream-function version keeps the central height to a more reasonable 
value. The treatment of the high off the west coast of the United States 
is similar. 



lli 



6. Conclusions 

The stratified model as proposed by ’Arnason appears, on the basis 
of the test cases presented here, to produce reasonable forecasts for 
levels at and beloi^ 500 mb. The stream-function version of the strati- 
fied model produces 500-mb forecasts for 2h and I48 hours which are in 
all cases very slightly better than those produced by the barotropic 
model. At 850 mb and 700 mb the model produces forecasts of system 
movement which are very realistic. The geostrophic version of the model 
suffers from the "blow-up" of low-latitude high pressure cells, as is 
common with most geostrophic models. 

Optimum values of the multiplying factor of the divergence term 
k]^ at 850, 700 and 500 mb havd been found to be 3.5? 2.75 and 0.9 
respectively. Although the number of cases tested was not large, the 
values determined here should be sufficiently accurate to provide a 
basis for the long-term testing necessary before any model becomes suit- 
able for operational use. 
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Figure U. Plot of vs. mean RMSE 
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Figure 850-mb analysis. 



20 January 1962 
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Figure 6. 850-mb analysis, ZL January 1?62 
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Figure 7» 8^0-mb analsrsis, 22 January 1962 
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Figure 8. 



2U-hour forecast for 8^0 mb from 20 January 1962, k^=0 



23 




Figure 2U-hour forecast for 8^0 mb from 20 January 1962, 
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Figure 10. 2ij.-hour forecast for 8$0 mb from 20 January 1962, 
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FigTore 11. US-hour forecast for 8^0 mb from 20 January 1962, k^=0 
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Figure 12. i;8-hour forecast for 850 mb from 20 January 1962, ^^=3 
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Figure 13. U3-hour forecast for 850 mb from 20 January 1962, k^=5 
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Figure lh» 850-mb divergence, 20 January 1962 
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Figure 15. 700-mb divergence, 20 January 1962 
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Figure l6. divergence, 20 January 1962 
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Figure 17. 1000-mb analysis, 6 1962 
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Figure l8. 8^0-mb analysis, 6 Ifey 1962 
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Figure 19. 2U-hour forecast for 1000 mb from 6 May 1962 



3it 




Figure 20. 2li-hour forecast for 8^0 mb from 6 l^ay 1962 
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Figure ZL. 1000-mb analysis, 7 I’fey 1962 
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Figvire 22. 850-mb analysis, 7 fey 19^2 
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Figure 23. U8-horar forecast for 1000 mb from 6 May 1962 
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Figure 2h» U8-hour forecast for 8^0 mb from 6 May 1962 



39 



I 




Figiire 2 $, 1000-mb analysis, 8 May 1962 
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Figure 26. 850-mb analysis, 8 May 1962 
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U November 1961 



- k2 



Figtire 27. 500-irib analysis. 




Figure 28. U8-hour forecast for 500 mb from I4. November I96I 
(stream function) 
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Figure 29. i;8-hour forecast for 500 mb from U November 1961 
(geostrophic) 
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Figure 30. 500-mb analysis, 



6 November 1961 
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